THE INTEGRALS IN GRADHTEYN AND RYZHIK. 
PART 5: SOME TRIGONOMETRIC INTEGRALS 



TEWODROS AMDEBERHAN, LUIS A. MEDINA, AND VICTOR H. MOLL 



Abstract. We present evaluations and provide proofs of definite integrals 
involving the function cos" x. These formulae are generalizations of 3.761.11 
and 3.822.1, among others, in the classical table of integrals by I. S. Gradshteyn 
and I. M. Ryzhik. 



1. Introduction 

The table of integrals [4] contains a large variety of evaluations of the type 

(1.1) ^~ I A{x)R{si-RX,cosx) dx 

J a 

where A is an algebraic function, R is rational and — oo < a < b < oo. We present a 
systematic discussion of two families of integrals of this type. This paper is part of 
a general program started in O |TOl [HI [1^ intended to provide proofs and context 
to the formulas in [4] . 

The first class considered here corresponds to the complete integrals 

.7r/2 

(1.2) c{n,p) := / xPcos'^xdx, 

Jo 

and 

1-11/2 

(1.3) s{n,p) / x^sin^xdx, 

Jq 

where n, p G N. In section [5] we present closed- form expressions for these integrals. 
These expressions involve the sums 

(1-4) E ' 



l<ki<k2<-<kj<n ^ ^ 3 

that are closely related to the multiple zeta values 

(1.5) C(«l,«2, • ■ • ,«s) = E 



The reader will find in Section 3.4 of |3| an introduction to these sums. 
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In general, one does not expect such elementary evaluations to extend to p ^ N. 
For example, the change of variables x = irt^ /2 produces 

(1.6) j x'^^^ cos xdx ^ V2tt J coa\^—jdt. 

The latter integral is evaluated in terms of the cosine Fresnel function 

(1.7) FresnelC[a;] := J cos f — j dt, 

which indeed is not an elementary function. 

The second class considered here presents generalizations of the formula 3.822.1 
in [4] stated as 

/■oo ^^g2n+l ^ 1 /F^ / 2n + 1 \ 1 

Jo ^^'"^^V2SU + ^- + JvlfTT' 



(1.8) 

The integral in (|1.8p can be transformed via t ^ to provide the evaluation of 
(1.9) / cos2"+ii2rft, 



fc=0 

„2 



Jo 

that is given as the case p = 2 in Theorem 13.21 

Section [3] contains analytic expressions for the generalizations 

/>oo 

(1.10) C„(p,6):=/ x-Pcos^''+\x + b)dx, 

Jo 

and 

/•OO 

(1.11) Sn{p,b):= x-Psm^''+\x + b)dx. 



The last section also contains some evaluations obtained by differentiation with 
respect to parameters. An illustrative example is 



oc />oo 



(1.12) / / log ^ log y + y)dxdy=ij + 2 log 2)^^ 
Jo Jo 

that is equivalent to 

poo poo 

(1.13) / / logx logy cos(x^ + y^) dxdy = — (7 + 21og2)7r^ 
Jo Jo 16 



A generalization of this evaluation appears as Example 13.31 

The method described in the present work gives impetus to a class of integrals 
that are closely related to the particular integral computations addressed in this 
paper. 

2. The first example 
In this section we present the evaluation in closed-form of the definite integrals 

(2.1) c{n,p) := / a:^cos"a;da;. 

"'0 

A special case of this appears as 3.822.1 in ,4j. 

The first step towards the evaluation of c{n,p) is to produce a recurrence. 



Theorem 2.1. The integral c{n,p) satisfies the recurrence 
(2.2) 

for n > 2, p > 2. 



(2.2) c{n,p) = !L^c{n-2,p)~ ^^P „ ^K {n,p-2), 
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Proof. The identity cos^ x — 1 ~ sin^ x yields 

(2.3) c{n,p) = c{n-2,p)- I{n,p) 

where 

r'r/2 



Now 



I{n,p) 



I{n,p) 



Tr/2 



COS X sm X ax 



d 

dx 



71-1 



cos" ^ x] dx 



1 

/ [x^ COS x + px^ ^ sin cos" ^ xdx 

~ ^ Jo 



2n 

c{n,p) ^ p 

71—1 71 — 1 



7r/2 



^ sin X 



cos" ^ X da;. 



Moreover 



7r/2 



x^ ^ sin X cos" X da; = 



7r/2 





p-1 



x'P ^4- I — — cos" X ) dx 

71 



dx 

c{n,p~2). 



□ 



Strategy: According to (|2.2p . the integral c{n,p) can be evaluated in terms of the 
initial values given in the table. The indices 77i and q have the same parity as n 
and p respectively and range over < m < n and < q < p. 



n modulo 2 


p modulo 2 


initial conditions 








c(77i,0) c{0,q) 


1 





c{m,0) c{l,q) 





1 


c{m,l) c{Q,q) 


1 


1 


c{m,l) c{l,q) 



We now evaluate the initial conditions c(7i, 0), c(7i, 1), c(0,p) and c(l,p). 

The expression for c(0,p). 

The computation of the identity 



(2.4) 

is immediate. 



c(0,p) 



1 /7r\P+i 



P+ 1 V2 



The expression for c(n,0). 

This is classical. The result appears as 3.621.3 and 3.621.4 in [4]. 



Theorem 2.2. (Wallis' formula and companion). Let ti e Nq. Then 

(2-5) ^(2-,0) = ^(f 

and 

22n 



(2.6) 



c(2ii+ 1,0) 



(2-+i)e:) 
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The shortest proof of Theorem 12.21 employs the representation 

(2.7) c(n, cos" xdx = (^^^, , 

that appears as 3.621.1 in [1]. Here B is the Euler's beta function defined by the 
integral 

(2.8) B{x,y)^ f t'^-^l - t)y^^ dt. 



The expression (|2.7p follows from the change of variables t = cos u. To express (12.5 
and (|2.6|) . in terms of the beta function, employ the standard relation 



(2.9) 5(.,,) = ^^, 

T{x + y) 

and the special values 

(2.10) r(n) ^ (n - 1)! and r(n + \) = ^^Ml 
that are valid for rt e N. 



The identity in Theorem 12. 2i in the case n is even, that is 

(2.11) ^(2n,0) = | cos^-ede = ^^[^^ 

is Wallis's formula and sometimes found in calculus books (see e.g. [5], page 492). 
To prove it, first write cos^ 9=1 — sit? 9 and use integration by parts to obtain the 
recursion 

2n — 1 

(2.12) c(2n,0) = c(2n-2,0). 

2n 

Then verify that the right side of (|2.1ip satisfies the same recurrence together with 
the initial value 7r/2 for n = 0. 

We now present a new proof of Wallis's formula (|2.1ip in the context of ratio- 
nal integrals. Extensions of the ideas in this proof have produced rational Landen 
transformations. The reader will find in [TJ [51 [51 [71 [S] details on these transforma- 
tions. 

Start with 

c(2n,0) - j cos2"0d0 = y ^ l + cosltj \ 

Now introduce ip = 29 and expand and simplify the result by observing that the 
odd powers of cosine integrate to zero. The inductive proof of (|2.1ip requires 

(2.13) c(2n,0) = 2- U]'^'^'^^^- 

Note that c(2n, 0) is uniquely determined by (|2.13p along with the initial value 
c(0, 0) = n/2. Thus (j^Tll) now follows from the identity 
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We now provide a mechanical proof of ()2.14p using the theory developed by Wilf 
and Zeilberger, which is explained in [131 [Mj; the sum in (I2.14|) is the example used 
in [2] (page 113) to illustrate their method. The command 

ci(binomial(n, 2i) binomial(2i, i)2^^*, 1, i, n, N) 

produces 

(2.15) /(n+1) = ^^/(n), 

71+1 

and one checks that 2^"(^^) satisfies this recursion. Note that (|2.12p and (|2.15p 
are equivalent under 



c(2n,0) = -^/(n). 



The proof is complete. 



Closed form expression for c{l,p). 
We now consider the evaluation of 

(2.16) c(l,p) / x^cosxdx. 

Jo 

The following evaluation appears as 3.761.11 in [3]. 

Theorem 2.3. Let p e N and (5odd,p be Kronecker's delta function at the odd 
integers. Then 

(2-17) c(l,p) = ^^(-1)^-^-^ (-) - (-l)«^5.,d,pp! 

where = [f J . 

Proof. Both sides of the equation ()2.17p satisfy the initial value problem 

TT- 2 

2 ■ 

Actually the recurrence (|2.18[) is obtained using integration by parts in (|2.16p . 
Iterating this recurrence yields the right hand side of (|2.17p . □ 



/7r\P 

(2.18) Up - p{p - l)up-2 ^ y-^j and lio = 1, ui 



Note 2.4. The result in Theorem 12.31 can be expressed in terms of the Taylor 
polynomial for cosx: 



(2.19) /,(a:) = (-l)S! -1 + E 



(-1) 



fe 



The formula (|2.17p can be restated 
(2.20) c(l,p) = J 



/p(7r/2), for p odd, 
/p(7r/2), for p even. 



2fe+l 



Closed form expression for c{n, 1): in fact, this would be the last initial condi- 
tion we require to execute the startegy outlined at the beginning of this section. 

Theorem 2.5. The integral c{n, 1) satisfies the recurrence 
(2.21) c(n,l) = !i^c(n-2,l)-i^. 
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Proof. The identity cos^ x = 1 — sin^ x yields 

(2.22) c(n, 1) = c(n - 2, 1) - J, 
where 

.7r/2 

(2.23) / a; sin^ a; cos""^ X dx. 



Integration by parts leads to 



1 



tt/2 



(2.24) J= / (sin a; + a; cos x) cos" ^xdx. 

n-1 Jo 

This produces (E^B)- □ 

The solution of (|2.2ip yields a closed-form formula for c{n, 1). 

Theorem 2.6. The integral c{n, 1) is given according to the parity of n, by 

/ 2 " o^*; \ 



(2.26) c(2n + l,l) = [ - -Y 

(2n+l)(2") I 2 



\ fc=l V fc / / 

for even indices. For odd indices, we have 

(- V O \ 

(2n+l)f;)\^2 A. 22fe (2fc + 1) y ■ 
To establish this result we solve a more general recurrence than (|2.2ip . 

Lemma 2.7. Let an, bn and rn be sequences with Un, bn ^ 0. Assume that z„ 
satisfies 

(2.27) a„z„ = bnZn-i + r„, n > 1 
with initial condition zq. Then 

(2.28) z„ = Mi:::^ f zo + T."T ^') 

aia2 ■■■an y ^ 6162 ■■■bk J 

Proof. Introduce the integrating factor dn with the property that fi„6„ — fi„_ia„_i. 
The recurrence (j2.27p becomes 

(2.29) Dn - Dn-i = dnrn, 
where Dn = (i„a„z„. Therefore, by telescoping, 

n 

(2.30) Dn^ Do + J2dkrk, 

fe=i 

with Dq ~ doaQZQ. To find the integrating factor, observe that 

(2.31) -/^ = 

"n-l On 

Thus 

(2.32) dn = 7^/--""-^ do. 

O1O2 ■■■bn 

Replacing in 1^:^ yields □ 
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Corollary 2.8. Let n G N and assume that z„ satisfies 
(2.33) 2nz„ ~ {2n — l)z„_i + r„, ti > 1, 

with the initial condition zq. Let Xj — 2^^ (^r') , then 



(2.34) 

Proof. Use a„ = 2n and fe„ = 2n — 1 in Lemma 12.71 □ 

We now apply Lemma l^TTI on the recurrence (|2.2ip . repeated here for convenience 
to the reader, 

n — 1 , , f 

c(n, 1) = c{n - 2, f) - — . 

n 

Observe that this recurrence splits naturally into even and odd branches. The 
value of c{2n, 1) is determined completely by c(0, 1), and c{2n + 1, 1) by c(l, 1). 
Hence, there is no computational interaction between c(2n, 1) and c{2n + 1, 1). Let 
Xn = c(2n, 1) so that x„ satisfies 

(2.35) 2nxn = (2n - l)xn-i - 
with the initial condition 

(2.36) a;o = c(0,l) = ^. 

Similarly, j/„ = c(2n + 1, 1), the odd component of c(n, 1), satisfies 

(2.37) {2n + l)yn = 2nyn-i - 

Zn + 1 

and the initial condition 

(2.38) yo = c(l,l) = |-l. 

The expressions for z„ in Lemma 12.71 yield the formulas for c(2n, 1) and also 
c{2n + 1,1) in Theorem 12.61 The proof is complete. 

Note 2.9. The finite sums in p.25p and (I2.26|) do not have closed-form, but it is 
a classical result that, in the limit, 

(2.39) y: 



and 



(2.40) J2 



2k\ 



22'=(2/c + l) 2' 



Note 2.10. Formula 3.821.3 in ^ gives formulas equivalent to ((2?25l) and (|2?26)) . 
respectively. 

Finally, we conclude this section by presenting ithe sought for closed form ex- 
pression for the integral c{n,p), for arbitrary n, p £ N. The recurrence (|2.2p . in the 
case of even indices n, becomes 

(2.41) 2nXnip) - (2n- l)X„_i(p) - ^fcllx„(p- 2) 

2n 
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where Xn{p) — c{2n,p). The initial value 
(2.42) Xoip) = 



1 



-TT- 



p+1 



(p+l)2P+i 

given in (|2.4p and the recurrence (|2.4ip show the existence of rational numbers 
0'n,p,p+i-2j such that 

(2.43) X„{p) - J2 an.p,p+i-2j7rP+i-2^', 

j=o 

with = [|J . The recurrence (|2.4ip is now expanded as 

(2.44) 2n^a„,p,p+i_2j^^'+i-'^' = (2n - 1) ^ a„_i,p,p+i_2j^^'+i-'^ 

j=o j=0 



an,p-2,p-l-2j"7r'' 



2n 

The fact that the coefficients Un.pj S Q allows us to match the corresponding 
powers of tt in (j2.44p . The highest order term is tt^^^. Only two of the sums contain 
this power, therefore 

(2.45) 2nan,p,p+i = (2n - l)a„_i,p,p+i. 

The initial condition 

1 



(2-46) «0'P^''+i = -(p+i)2P+i 

comes from (j2.42p . The solution to the initial value problem (|2.451 12.46P is then 
found using Corollarv l2.8l (here r„ = 0), namely that 

a„,p,p+i = 2)22n+p+i ■ 
The coefficient of the next highest power n^^^, in (|2.44p . yields the recurrence 

pip — 1 ) 

(2.48) 2na„_pp_i = (2n — l)a„_i_pp_i — — an,p-2,p-i- 

Observe that the last term in this relation is given by (I2.47p . Moreover, (|2.42p 
shows that ao.p,p_i = 0. The solution to (|2.48p . following Corollarv l2.8[ is 

p^2jij n ^ 

(^■^^) an,p,p-l = ~ 22"+P+l ^ l2' 

fcl=l 1 

The next power of tt in (|2.44p produces 

p{p-l){p-2) f2n 



(2.50) 2nan,p,p-3 = (2n - l)a„_i,p^p_3 + 



j22n+p 



f2n\ ^ 1_ 

[nj ^ 

^ ^ fci=i 1 



with ao,p.p-3 = 0. One more use of CoroUarv 12.81 vields 

(2-51) an,p,p-3 = 22"+p+i(p-3)! ^ ^ k^' 

This procedure can be repeated until all descending powers of tt have been ex- 
hausted, hence a complete closed form for the integrals c{n,p) will be made possible. 



TRIGONOMETRIC INTEGRALS 



9 



Theorem 2.11. Let n, p e N and let = [|J . Then the even branches Xn{p) 
c{2n,p) of the integral 

(2.52) c(n,p) / xPcos^xdx 



JQ 

are given by 

(2.53) Xnip) = ^ a„^p^p+i„2j7r''+^"^-' + (5odd,p ■ a^,p, 

and the value of an.p,p+i-2j for P > 2 and < j < is given by 

an,p,p+l-2j - 2„+p+i {p+l- ^ klkl • • • fcf 

•'^ l<fci<fe2<---<fcj<n 12 J 

and 

(2 54) a* -tltf>[ y 1 Y 

l<fci<fc2<---<fcp<n 1 2 p j = l Vj/ 

Similarly, for the odd branches Y{n^p) = c{2n + l.p) we have 

(2.55) r„(p) = 6„,p,p-2j7rP-2j + j^^^ p . 

i=o 

with 

(„l)jp!22»+2j-P 1 



(2n + !)(-) (p - 2,)! „<,^<,^..<,^.<„ (2fci + l)^(2fe + 1)^ • • • (2/c, + If ' 



and 

, (-l)^^p!22" ^ ^ ^ 

(2n + o<..<.^..<..<„ + l)^(2fc2 + 1)^ • • • (2/cp + 1)2 A. 22.(2, + 1) ' 

3. Some examples on the halfline 
In this section we provide an analytic expression for 

(3.1) C„(p,6)=/ x~Pcos^''+^{x + b)dx, 

Jo 

and 

/"OO 

(3.2) Sn{p,b)= x-Psm'^''+\x + b)dx. 

Jo 

In the table [3] the evaluation of the special case P = \ and 6 = 0: 
cos2"+ix ^ 1 /F^/2n + l\ 1 



' 22" V 2 ^ V'^ + ^' + V %/2fcTT' 

and 

psin2"+ix^ 1 /FJ^, 2n + l \ 1 







X 22" V 2-^^ ^ V'^ + ^ + 1/ V2fc + l' 



as 3.822.2 and 3.821.14. 
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Theorem 3.1. Let < p < 1 and n e No := N U {0}. Then 

for\ f°° -P 2„+l . r(l -p) . fTTp\^ Cn-k) 

(3.5) / X ^ cos ^ xax — sm — > — — 

^ ' Jo 22" V 2 / ^ {2k + iy-P 



and 



,na\ f°° -p ■ 2«+i , r(l-p) /7rp\ -A ^ Cn-k) 

(3.6) / X ^ sm ^ xdx= — — cos — > —1 7-- -r- ■ 

^ ' Jo 22" \2J^^ ' (2A; + l)i-P 

Proof. The identity 2 cos a; = + e~" and the binomial theorem yield 
(3.7) 

X'PcOS^'^+'xdx = 2-2"-l J2 (^''^ ^) J^^ (e»(2"+l-2fe)- + g-.(2„+l-2/c)^ 



k=0 



Recall the Heaviside step function defined by H{x) = 1, if a; > and H{x) — 
otherwise. Then, each of the integrals in (|3.7p is evaluated using the Fourier 
transform 



(3.; 



H{x)x~Pe-"^'' dx = ^jL|^exp(-z7r(l -p)sign(cj)/2). 

I 1^1 



□ 



Corollary 3.2. Let p > 1 be real and n <E No. Then 

n f2n+l\ 



(3.9) 

and 
(3.10) 



r COS-+1 x^ dx^^vi ^-±l\ cos ( ^\ y , 
Jo 22" \ p ) \2pJ f^^{2k + iy/p' 



(2n+l\ 
: \n-k) 

(2fc+ l)i/p- 



Proof. The change of variables x 1-^ x^^^-^ in the results of Theorem 13 . 1 1 gives the 
result. □ 

The last result described here is a further generalization of Theorem 13.11 
Theorem 3.3. Assume 6eM, 0<p<l and n e No. Define 

/•OO 

(3.11) C„(p,6)= / x^P cos^''+\x + b) dx 



and 

/>OC 

(3.12) Sn{p,b)^ x-Psm'^''+\x + b)dx. 

Jo 

Then 

f^^-i^ r( n r(l-rt ^^2» + l^ sin(f -(2fc + l)b) 



and 



,3., ,,,.,.£0^ -'I;'-:"-' 
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Proof. Denote the left-hand side of p.l3|) and (|3.14|) by fn{b) and gn{b) respectively. 
Differentiation with respect to the parameter b yields 

(3.15) ^-(-ir(2n+i)/„ = {2n+i)j2i-^yC')m 

^ + (-ir(2n+l)g„ ^ -{2n+l)J2{-^y(fj9Ab)- 

Considering b and p fixed, wc now show that the right-hand side of p.l3p and p.l4p 
satisfy the system p.lSp with the same initial conditions. This will establish the 
result. 

In the case of the right-hand side of (|3.13p . it is required to check the identity 

"2n-t- 1\ sin(fp- (2fc+ 1)6) 



{2k + iy-p 

2j + l\ sin(fp- (2fc-h 1)6) 



(2n-fi)y(-i)-' " 2-2^y, , , , 

To verify this we compare the coefficients of the transcendental terms 

sm{^p - {2k + l)b) 

{2k + iy-p ■ 

It turns out that this question is equivalent to validating the identity 

(3.16) (-1)'^2-" i2k + 1) = (2n + 1) ^(-1)^2-^ (^) (^^ + ^ 

j—k 

To this end, we employ the WZ-technology as explained in [14, . This method 
produces the recurrence 

(3.17) 2{n + k + l){n + 1- k)u{n + 1 - fc) - (n + l)(27i -f- 3)u(n, k) = 0. 

To prove p.l6p simply check that both sides of p.l6p satisfy the recurrence p.l7p 
as well as the initial condition u(0,0) = 1. 
The identities 

/•OC 

(3.18) / x^P cos{x + b) dx = -r(l -p)sin(6- 2f^) 

/ x~P sin{x ^ b) dx = r(l - p) cos(6 - ^), 
Jo 

which are special cases of 

(3.19) / x-P cos{ax + b) dx = -aP-^r{l ~ p)sm{b ~ l^) 
Jo 

x~P siniax + b) dx = aP-^T{l^p)cos{b-If), 



show that the corresponding initial values in (I3.13P (respectively [3TT3]) match. The 
evaluations p.l9p appear as 3.764.1 and 3.764.2 in [1]. To establish (I3.18P expand 
cos(x-|-6) as cos a; cos 6— sin a; sin 6, use the change of variables a; i— > x^, and Theorem 
O □ 

We now discuss some definite integrals that follow from Theorem 
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Example 3.1. DifFcrcntiating ()3.13p with respect to p and setting p = ^ and 6 = 
gives, after the change of variables a; i— > a;^, 

(3.20) 

A f2n + 1\ 1 



poo r— 

/ logxcos2"+ia:2dx = -^^(7r + 27 + 41og2)^ 

" k=0 



^ f2n + 1\ \og{2k + 1) 



k J y/W- 



22"+2^V"-^/ V4fcT2 ' 
where we have used the value r'(l/2) = — \/7r(7 + 2 log 2). 

Example 3.2. Assume < p, q < 1. Multiplying p.l3p by and integrating 
over the half-line yields (after replacing b by y) 

r r ^''y^-y) ,A - -r,i-p)r,i-,)c„s^'<'' + ') 

Jo JO ^ y 

n 

X E 

In particular, for n = 0, 



2 

2n + l\ (2fc + l)P+9-2 



n-k J 22^^ 



OO /'OO 



The derivative g^g^ at p = = ^ produces the evaluation 

(3.22) / / log ^ log y ^^Qg/ +y)dxdy^h + 2 log 2)7r^ 



that we promised in the Introduction. 

Example 3.3. Iterating the method described in the previous example yields 

f ( II ||2N TTi (-l)'^"7r"'/^ /ReV'n if n is even, 

Jk" 1 ^ Im V'n if IS odd, 



with 



(3.23) A„ = ^„ = + 2 log 2 + ^) " e-/^ 

Here WxW^ ~ x\ + ■ ■ ■ + x\ and 7 is Euler's constant. For instance, for n = 3 we 
have 

log X logy log zcos(x^+y^+z^) dec dydz = — — (-16$^+12^^7r+6^7r^- 



/O JQ JQ 

where ^ = 7 + 2 log 2. 
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